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Abstract 

In this paper, we consider the initial and boundary value problem of a simplified nematic 
liquid crystal flow in dimension three and construct two examples of finite time singularity. The 
first example is constructed within the class of axisymmetric solutions, while the second example 
is constructed for any generic initial data (uo,do) that has sufficiently small energy, and do has 
a nontrivial topology. 


1 Introduction 

Let n C M” (n = 2, 3) be a bounded, smooth domain, and 0 < T < +oo. In this paper, we will 
consider a simplified version of the hydrodynamic flow of nematic liquid crystals on 0 x (0, T) given 
by 

ut + u ■ Vu — pAu + VP = —AV • (Vd © Vd — ^ j 

< V-u = 0, (1.1) 

+ u ■ Vd = 7(Ad + |Vdpd), 

where u{x,t) : 0 x ( 0 ,r) ^ M” is the velocity field of the underlying incompressible fluid, d(x, t) : 
n X (0,T) ^ := {u G : |u| = l} represents the (averaged) orientation field of nematic 

liquid crystal molecules, P(x, t) : U x (0, T) —>■ M is the pressure function, x G U, V- denotes the 
divergence operator on M”, Vd 0 Vd = ((^i represents the stress tensor 

induced by the orientation field d, and is the identity matrix of order n. 

The parameters fi, A and 7 are positive constants representing the fluid viscosity, the competition 
between kinetic energy and potential energy, and the macroscopic elastic relaxation time for the 
molecular orientation field respectively. 

The system mil, first proposed by Lin [2S], is a simplified version of the general Ericksen-Leslie 
system modeling the hydrodynamic flow of nematic liquid crystal materials proposed by Ericksen 
m and Leslie [21] during the period between 1958 and 1968. The system mi is a macroscopic 
continuum description of the time evolution of the material under the influence of both the fluid 
velocity field and the macroscopic description of the microscopic orientation configurations of rod¬ 
like liquid crystals. The interested readers can refer to m, m, m, and Lin-Liu m for more 
details. 

The system dni) is a strongly coupling system between the incompressible Naiver-Stokes equa¬ 
tion and the heat flow of harmonic maps into §^, and relates to several important equations: 
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1 . When A = 0 or d = eo G the system 01,2 reduces to the incompressible Navier-Stokes 
equation (or NSE) that has been extensively studied for decays (see Lions [35], Temam [39]L 
Although the existence of global weak solutions to the initial value problem to NSE has been 
established by Leray [23] in 1930’s and Hopf [20] in 1950’s, it is a long outstanding open 
problem whether NSE admits a global smooth solution for a smooth initial data in dimension 
n = 3. It is also an open problem whether Leray-Hopf weak solutions to NSE are unique in 
dimension n = 3. 

2. When u = 0 and 7 = 1, the equation ms reduces to the heat flow of harmonic maps into 
For dimension n = 2, the existence of a unique global weak solution, which has at most finitely 
many singular points, has been proved by Struwe [38] and Chang [T]. In higher dimensions, 
the existence of a global, partially regular weak solution has also been obtained by Chen- 
Struwe [1] and Chen-Lin [5]. Examples of finite time singularities have been constructed by 
Coron-Ghidaglia [7] and Chen-Ding [ 6 ] for n > 3. In an important work [3], Chang-Ding-Ye 
constructed examples of finite time singularities when n = 2 (see Grotowaski mm for some 
generalizations to n > 3) by studying the equation 

iPr sin( 2 c?) / s / ^ , X . s 

V>t = ^rr + - -(Gt) G (0,1) X (0,+oo). (1.2) 

r 2r^ 

The interested readers can refer to Lin-Wang m and references therein for more details. 

3. Another important case we want to mention is 7 = 0. The system m for (u, d) is closely 
related to the MHD system for (tt, ip) provided we identify ijj = X/ x d. There have been many 
interesting works on global small solutions to the MHD system recently. See, for example, 
Lin-Zhang [32] and Lin-Zhang [33] for n = 3 and Lin-Zhang-Xu [3l] for n = 2. 

The system (II.Ih has attracted a lot of interests and generated many interesting research works 
recently. Here we would like to mention a few of previous results. In dimensions two, Lin-Lin-Wang 
|26j have proved the existence of global Leray-Hopf type weak solutions to (II. ip with initial and 
boundary conditions, which is smooth away from finitely many possible singular times (see Hong 
m for (II.Ih in D = R^, Hong-Xin [18] and Xu-Zhang [41] for other related works). Lin-Wang [29] 
have also proved the uniqueness for such weak solutions. It remains a very challenging open problem 
to establish the existence of global Leray-Hopf type weak solutions and partial regularity of suitable 
weak solutions to (jl.ip in dimension three. Very recently, Lin-Wang [28] have proved the existence 
of global weak solutions in dimension three under the assumption that do{x) G for a.e. x G Cl 
by developing some new compactness arguments, here is the upper hemisphere. When H = R^, 
the local well-posedness of m was obtained for initial data {uq, do), where {uq, Vdo) G Lj)[Q^(R^), 
the space of uniformly locally L^-integrable functions, has small norms, by Hineman-Wang m- 
While the global well-posedness of (II.ip was obtained by Wang [iO] for (uo,do) G BMO ^ x BMO 
with small norms. A BKM type blow-up criterion was obtained for local strong solutions to m 
by Huang-Wang m (see also Hong-Li-Xin [l9]). More references can be found in the survey paper 
by Lin-Wang [30] . 

It is a very interesting question whether the short time smooth solutions to the nematic liquid 
crystal flow m develop singularities in finite time. It is not hard to verify that in dimension two, 
if d is the heat flow of harmonic maps with finite time singularity, constructed by [3] , and if we set 
tt = 0, then (u, d) is also a solution of (11.11) which has a finite time singularity. However, it will 
be more desirable to construct an example in which the fluid velocity field u is non-trivial. In this 
paper, building upon the construction by [3] on the heat flow of harmonic maps, we are able to 
construct in dimension three the first example of solutions to ( 11 . 11 ) with finite time singularity in 
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which both the fluid field u and the director field d are non-trivial. More precisely, we consider 
in the class of axisymmetric solutions without swirls and show that for a suitably chosen domain 
Q and initial-bounday data {uq, do), the short time smooth solution {u, d) to (jl.ll) develops a finite 
time singularity. 

Let C M”" denote the unit ball centered at 0. Now we state our first result. 

Theorem 1.1 Let Vt = x [0,1]. There exists (po G C'°°([0,1]), with ¥^ 0 ( 0 ) = 0 and |¥5o(l)| > tT; 
such that if 

uo(x) = {x,y,-2z), 

and 

= ( I sinv?o(Vx^+ y2), — siny)o(Va;^+ ^2^,cosv?o(\/a:^+ ?/^)), 

for X = (x, y, z) G Ll, then the short time smooth solution {u, d, P) to the system lll.l]) in Q, under 
the initial and boundary condition: 



(u(x, 

0),d(x,0)) = 

{uo{x),do{x)), X 

G LI, 

(1.3) 


u(x 

,t) = 'Uo(x) 

X G dQ, t > 0, 



< 

d(x 

,t) = do{x) 

X G dBf X [0,1], 

o' 

A 

(1.4) 


.mi- 

-,t) = o 

X G X {0,1}, 

t > 0 . 


must blow up at time Tq for 

some 

0 < To = ro((/?o) < +00. 




Since is axisymmetric and (uo,do) is axisymmetric without swirls, the uniqueness of short 
time smooth solution [u, d) of (|l.ip implies that it is axisymmetric without swirls. By converting 
the system (jl.ip into the form of being axisymmetric without swirls, the proof of Theorem 11.11 
utilizes two interesting observations: 

1) The velocity u is a static solution to the Navier-Stokes equation. In fact, u{x,y,z,t) = 

(x, y, —2z) is spatial gradient of the quadratic harmonic polynomial h{x, y, z) = + — 

2) The angle function (p{r, t), associated with the orientation field d(x, t), solves a drifted version 
of the equation (|1.2p : 

iPr sin(2(^) , s . X , N 

Pt + rpr = Prr + ^ -2r^’ ^ ^ 

After a suitable re-parameterization to handle the contribution from the drifting term rprii", t) 
in we can modify the construction of [S] to build a subsolution to (|1.5I) that blows up at 
finite time. This, combined with the comparison principle, yields the finite time singularity 
of (fL5l) . 

In particular, u is smooth in this example. This is consistent with known results of the Navier- 
Stokes equation: any local axisymmetric solution to the Navier-Stokes equation, without swirls, is 
globally smooth (see, for example, Leonardi-Malek-Necas-Pokorny [22]). Therefore, the finite time 
singularity arises essentially from the orientation field d. 

Now we would like to make a few comments related to Theorem ll.il 
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Remark 1.2 a) By modifying the argument of Chang-Ding [2], we will show, in Theorem 2.5 
below, that (II.5|) admits a global smooth solution if the initial-boundary value ipo G 1]) 

satisfies (^o(O) = 0 and ||v^o||(^qo ij) ^ 

b) It is a natural question to ask whether the stress tensor S := Vd © Vd — -|Vd|^l 3 blows up in 
finite time, if (u, d, P) is the solution constructed by Theorem 11.11 The calculations in Section 5 
seem to suggest that S may not blow up. 

c) Since uo(x) ■ u ^ 0 on d^l, here v denotes the outward unit normal of 911, it is unclear whether 

the solution {u, d, P) constructed by Theorem 11.11 enjoys the energy dissipation inequality for 0 < 
t < Tq: ^ 

4 [ + |Vd|2)(t) + 2 [ [ (/r|Vu|2 + -|Ad + iVdpdH < 0. (1.6) 

at Jn Jo Jn 7 

See Section 5 for more details. 

It is the question from Remark 1.2 c) that motivates us to construct another example in which 
the solution (u, d, P) to (jl.ll) develops finite time singularity, and satisfies the energy dissipation 
inequality (jl.6jl . 

In order to state it, we need some notations. Denote the north pole by e = (0,0,1) G Set 


:= [v G 


V • u = 0 


}. 


and 


C~(R)*,§2) := |d G (7°°(R?,S2) I d = eon SR?}. 


For continuous maps f,g £ C'(Rf,S^), with f = g on dBf, we say that / is homotopic to g 
relative to dBf if there exists a continuous map G C[Bf x [0,1],S^) such that 

(i) = /(•) = g(-) on dB?, for all 0 < t < 1; and 

(ii) $(-,0) = /(•) and $(•, 1) = g{-) in Bf. 

Now we have 

Theorem 1.3 There exists eo > 0 such that if uq G M^) and do G satisfies 

that do is not homotopic to the constant map e : Bf —?■ relative to dBf, and 


E{uo,do):=l [ (|uoP + |VdoP) <e§. 


2 JbI 

Then the short time smooth solution {u, d, P) : Bf x [0, T) 
flow lll.l]) . under the initial-houndary condition 


(1.7) 


X to the nematie liquid crystal 


(w,c^)|i=o = (^0,4), in ^1, 

{u,d)\g ^3 = (0,e), 0 < f < T, 


( 1 . 8 ) 


must blow up before time T = 1. 

The following remark indicates that there are ample examples of {uo,do) G {Bf,R^) X 

C^{Bf,§'^) satisfying the conditions of Theorem 11.31 

Remark 1.4 a) Let H{z,w) = (|zp — \w\‘^,2zw) : = {{z,w) G C x C : \z\'^ + \w\^ = l} —)• C 

M X C be the Hopf map. Let Dx{x) = Ax : ^ be the dilation map for A > 0, II : —)• 
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be the stereographic projection map from e, and 'I'a = II ^ o D\ o II : ^ Then direct 

calcnlations imply that the Dirichlet energy oi H o'^\ : ^ satisfies 

lim [ \V{H da = 0. 

is® 

Moreover, it is easy to see that H o is not homotopic to the constant map e : ^ S^. Let 

$ G §3) such that {e} is a diffeomorphism and $ = e on dBf. Now we can 

check that for any u G since 

lim E{X~^u, H o Ta o $) = 0, 

A—^ oo 

we can find a sufficiently large Aq > 0 depending on u, H, and such that 

(no, do) := (Aq ^n, H o Taq o ^) : Bf ^ x 

satisfies the condition (|1.7I) of Theorem 11.31 and do is not homotopic to the constant map e relative 
to dBf. 

b) The initial data do constructed in a) has previously been used by Ding-Wang [9] in the con¬ 
struction of finite time singularity of the Landau-Lifshitz-Gilbert system modeling the continuum 
theory of ferromagnetism. 

Since precise values of the parameters A and 7 in (|l.ip don’t play a role in this paper, for simplicity 
we assume henceforth that 

A = 7 = 1. 

The paper is organized as follows. In Section 2, we will derive the axisymmetric form of dni), 
without swirls. In Section 3, we will sketch the proof of local smooth axisymmetric solutions of ()2.2p 
without swirls. In Section 4, we will prove the existence of global smooth axisymmetric solutions of 
(ESI) without swirls, when the initial data (po satisfies </?o( 0 ) = 0 and |(/?o(?’)| < for all 0 < r < 1 . 
In Section 5, we will present in details the example of solutions to (12.2p with finite time singularity 
for suitably chosen initial data (po G 1]), with (/Jo(0) = 0 and |<y7o(l)| > tt. In Section 6, we 

will outline the proof of Theorem 11.31 

2 Axisymmetric form of (ll.ll) without swirls 

In this section, we will derive the axisymmetric form of (jl.ip without swirls. We would like to 
mention that Dong-Lei [ 8 ] have constructed a global smooth axisymmetric solution of dLH) in 
dimension two. 

Let (r,9,z) denote the cylindrical coordinates of M^, and set 

e'’= (cos d, sind, 0 ), = (—sin 0 , cos d, 0 ), e^ = ( 0 , 0 , 1 ) 

as the canonical orthonormal base of in the cylindrical coordinates. For a G [0,27r], let Ra G 

SO(3) denote the rotation map of angle a with respect to the z-axis. Recall that a vector field 
— 7 - is axisymmetric if 

R~^ o V o Ra = V, V a G [0, 27r]. 

Hence any axisymmetric vector field v can be written as 

v{r, 9, z) = v^{r, z)e’' + v^(r, z)e^ -|- v^(r, z)e^. (2.1) 
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If, in addition, = 0, we say v is axisymmetric without swirls. 

A solution (u, d, P) of the nematic liquid crystal flow equation (jl.lj] is said to be axisymmetric 
without swirls, if 

u{r, 6, z, t) = u’'{r, z, t)e^ + u^(r, z, t)e^, 

< d{r, 6,z,t) = sin ip{r, z, t)e^ + cos z, t)e ^, 

^P{r,e,z,t) = P{r,z,t). 

A domain C is axisymmetric if it is invariant under a rotation map Ra for any a G [0,27r]. 
Now we have 


Lemma 2.1 For any axisymmetric domain Q C if {u,d,P) is an axisymmetric without swirl 
solution of the system ([HP in ^ X M_|_, then {u^ ip, P) solves 


where 


Proof. 


' Du'' 

15t 




pAu'' + ■^u'' + Pr = — (A(/5 

Du^ X 3 n /X 
- pAu^ + Pz = -{Aif 

-{ru'')r + iu^)z = 0, 

T _ 


sin(299) ^ 
2r2 

sin(2(/9) X 
2r2 


_ A = 

Dt 2r2 


' p_ 

< Dt 
A 


dt + u''dr + u^dz, 

dl + -dr + dz- 

r 


Let’s first show (j2.2h -i. By (II.ip ^ and the definition of u, we have 


( 2 . 2 ) 


0 = 9. + dy + dz (u^) 

_ dx{u''x) + dy{u''y) u''{xd^r + ydyv) ^ ^ ^^ 3 ,^ 

— 9 \ Uz\U ) 


2u^ 

r 


r\ r ^ 

=-h drU +Uz -= 

^ r 


{ru'')r 


+ ul. 


For (|2.2p .i. since d can be written as 


d = (cos 6 sin ip, sin 6 sin ip, cos ip), 


direct calculations imply that 

/ 

dt = (cos 9 cos ip, sin 6 cos ip, — sin ip)ipt, 
dr = (cos 9 cos ip, sin 9 cos ip, — sin ip)ipr, 
dg = (—sin 0 sin(/9, cos 0 sin(/9,0), 
dz = (cos 9 cos ip, sin 9 cos ip, — sin ip)ipz , 

and 

drr = {fcosipiprr — sin (cOS 0, sin 0), —{simpiprr + COS(/?(/9^)) , 

< dgg = (— cos 9 sin ip, — sin 9 sin ip, 0), 

^dzz = {{cosipipzz - smipipl){cos9,sm9),-{simpipzz + cosipipD) . 
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Hence we have 


IV7 2 2 I ^ ■ 2 I 2 

I Va| — H—2 Sin + <y9^, 

and 

Ad + IVdpd = ((/?rr + ^zz H—-^^~ ‘ 

By the definition of u, we also have 

u • Vd = {vi'dr + v?dz)d = {vJ'(pr + u^Pz) (cos 6 cos tp, sin 6 cos p, — sin p). 

Putting these identities into (ll.ip -?. we obtain 

„ o / 1 sin(2w)x 

Pt “b 11 Pr “b 11 Pz — \Prr “b Pzz “b 2^ ' 

This yields the equation (|2.2L . Now we want to derive the momentum equations (I2.2p i o. Since 

Ae' = -le', (2,3) 

it follows from 3ilici (|^.3p 

- pAu^ + ^u!' + Pr)e^ + - pAu^ + Pz)e^ = —Ad • Vd. (2-4) 

^ Dt Dt ' 

From (|2.3p . we also have 

Ad = (A- 2 ) (sin p)^^ + A(cos p)e^. 

Since 


Vd = V(sin(/9)e^ + V(cos(/9)e^ + sin(/9Ve^, 


and 


(Ve^e^) = (Ve^e3) = 0 , 

we have 

—Ad • Vd = — (a- 2 )(sinv 5 )V(sin<^) — A (cos 99 ) V (cos (^) 

= (A(cos p) sin (/? — (A- 2 ) (sin(^) cos p) {pr^^ + Pz^^) (2-5) 

, 1 sin(2(/5). , „ ^ o. 

— \Prr “b Pzz “b ~Vr 2^ (V^f® ~b Pz^ )■ 

Putting (j2.5|) into (j2.4l) yields (ll.ip i 9. □ 

In order to construct a solution of (12.2p with finite time singularity, we further consider the 
domain to be the round cylinder Q = x [ 0 , 1 ] and an axisymmetric, without swirl, solution 
(u, P, d) in the special form: 

n(r, 9, z, t) := v{r, t)e'" + w{z,t)e^, 

< d{r,6,z,t) := sin p{r,t)e'^ + cos p{r,t)e^, 

^P{r,e,z,t) := Q{r,t) +R{z,t). 
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Then (|2.2I) becomes 


( 1 \ ^ / 'fr sin(2(/3) ^ 

Vt + VVr - fJ.[Vrr + “ “ + Qr = (V^rr + — - r G [0, IJ, 

Wt + WWz — VWzz + = 0, z G [0,1], 

^{rv)r + Wz = ^, (r, z) G [0,1] X [0,1], 

(Pr Sin(2<y9) 

Pt + VPr = Vrr + — -2r^’ '^^[0,1]. 

The initial condition for (n, w, (p) reduces to 

v\t=o = r, 0 < r < 1, 

< tt)|t=o = —2z, 0 < z < 1, 

= ^o{r), 0 < r < 1, 

for some po G C°°([0,1]), with (/?o(0) = 0. The boundary condition is 


II 

0 

0 

II 

II 

0 

\v{l,t) = 1, 

[u;(l,t) = -2, 

\p{l,t) = (^o(l) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


Remark 2.2 It is easy to verify that the initial and boundary condition (I2.7l) - (|2.8p for {v,w,(p) is 
equivalent to the initial condition ()1.3p and the boundary condition (11.411 for {u,d). 

Now, we are ready to state our main results on the system (|2.6p . under the initial and boundary 
condition ()2.7p and (12. 8 p . 

The first result asserts that for any short time smooth solution (v, w, p, Q, R) to (|2.6p along 
with the initial-boundary condition (j2.7h and (j 2 . 8 h . {v,w) is static and Q,R are also static (up to 
a time-dependent constant). In fact, we have 

Lemma 2.3 For 0 < T < -|-oo, suppose that v,w £ C'°°([0,1] x [0,T)) satisfies 

^{rv)r + Wz = 0, (r, z) G [0,1] X [0,1], (2.9) 

and 

fn( 0 ,t) = 0 fu;( 0 ,t )=0 

|n(l,t) = l, |u;(l,t) = - 2 . 

Then v{r,t) = r for any (r, t) G [0,1] x [0,T), and w{z,t) = —2z for any {z,t) G [0,1] x [0,T). 
Proof. Differentiating (12.911 with respect to z yields 

Wzz{z,t) = 0 

so that w{z,t) = ai{t)z + a 2 {t) for some functions afit) and a 2 {t). Since w{0,t) = 0 and w{l,t) = 
—2, we see that 02 ( 1 ) = 0 and oi(t) = —2. Thus w{z,t) = —2z. 

Similarly, differentiating (12.9p with respect to r yields 

(“Mr-)^(Di) = 0, 

this implies that rv{r,t) = bi{t)r‘^ + h 2 {t) for some functions hi{t) and h 2 {t). Since n(0,t) = 0 and 
v{l,t) = 1, we see that b 2 {t) = 0 and 61 (t) = 1. Thus v{r,t) = r. The proof is complete. □ 


Next we have 








Lemma 2.4 For 0 <T < +oo, {v,w,(p,Q, R) E C'°°([0,1] x [0,T)) solves 112.6[) . under |g.7| j and 
(KM) iffv{r,t) = r, w{z,t) = -2z, 


R{z, t) = —+ ci{t) 

s r , ^Pr sin(2f/?). 

Q{r,t) = - {'Prr + ^ - ^^^)iprdr - — + C 2 it), 


for some ci,C 2 E C'°°([0,T)), and 

' ifj. sin(2(/5) 

(ft + npr = iPrr-\ -0 < r < 1, 

' (/?(r,0) = p>o{r), 0 < r < 1, 

= 0, = v9o(l), t > 0. 

Proof. Applying lemma [23] to (12.61) 9 yields that 

Rz + 4z = 0, 


( 2 . 11 ) 


( 2 . 12 ) 


this, by integration, implies 


R{z,t) = -2z‘^ + ci(t), 

for some function ci E C'°°([0,T)). Apply lemma 1^31 to (|2.6p i /t. we obtain 

_ . (Pr sin(2v;). 

Qr — \^rr H“ ^ ^^2 j 

ifr sin(2(/:?) 

“1“ “1 9 • 

T Zt^ 

Integrating (j2.13p yields (j2.1ip i . while f|2.1ip 9 follows from (j2.14p . (j2.7P ' and f|2.8p . 


(2.13) 

(2.14) 
□ 


It is readily seen that whether (j2.6l) . under (|2.7n and (|2.8I) . admits a global smooth solution 
is equivalent to whether (|2.12p admits a global smooth solution. For this, we are able to extend 
Chang-Ding [2] and Chang-Ding-Ye [3] on the heat flow of harmonic maps and obtain the following 
two results. 


Theorem 2.5 Suppose cpo £ ^““([0,1]) satisfies </?o(0) = 0 and |75o(’')| < for all r E [0,1]. Then 
there is a unique, global smooth solution {v,w,ip) to p2.6p . under (12.7p and (j2.8p . 

Theorem 2.6 There exists fio E C°°([0, Ij), with 4>o{0) = 0 and |(/9o(l)| > tT; such that the short 
time smooth solution {v,w,(p) to (12.6p . under (|2.7I) and (12.8p . must blow up at Tq for some 0 < 
To = To{(j)o) < + 00 . More precisely, ipr{0,t) —oo as t ^ Tffi. 

Remark 2.7 It is readily seen that Theorem 11.11 follows directly from Theorem 12.61 Since p>{r,t) 
blows up at r = 0 and t = Tq, we see that d{r,9,z,t) blows up at r = 0, z E [0,1], and t = Tq. 
Thus the singular set of d at the first singular time Tq includes the line segment {0} x [0,1] C 11. 

The proofs of Theorem 2.5 and Theorem 2.6 will be given in Section 4 and Section 5 respectively. 
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3 Existence of short time smooth solutions 


In this section we will establish the existence of short time smooth solutions to the system (j2.6p . 
under the initial and boundary conditions (j2.7ll and (12.8p . By lemma [231 it suffices to prove the 
local solvability of the drifted-harmonic map equation (I2.12h . 


Theorem 3.1 (local existence) For any ifo € C'°°([0,1]), with (^o(O) = 0, there exist to > 0 and a 
unique smooth solution ^p{r,t) on [0,1] x [0,to) lo the equation 112. 


Proof. First observe that if ip solves (12.121) . then 


d{x,y,t) := ( siny?(y/x^ + ^ ,sm(p{y/x‘^ + ,t) ^ - ,cos(p{y/x‘^ + y'^,t)), 


^x^ + y^ ' ^ ^ ^ -y/x^ + 

for (x, y) € and t > 0, solves 

dt + (x, y) • Vd = Ad + |Vd|^d, (x, y) G B\, t > 0, 


(3.1) 


with the initial condition 

d\^^o = (smipo{y/x^ + y'^) —^ -- ,sinyo(y/a:^ + y^) ^ cos ipo{Vx^ + y"^)) , (3.2) 

^ y x^ + y^ y x^ + y^ ^ 

and the boundary condition 

^IdBf = (sinyo(l)a;,sinyo(l)2/,cosyo(l))- (3.3) 

By the standard contraction mapping theorem, we can prove that there exist to > 0 and a unique 
smooth solution d G C°°{Bl x [0,to))§^) to (13.Ih . along with (13.2p and (13.3p (see, e.g., [26] Section 
3 for a detailed proof). 

Now we need to argue that such a solution d(x,y,t) to (j3.ip - (l3.2h - (l3.3p is axisymmetric, i.e., 

d(x, y, t){:= d{r, 9, t)) = (sin ip{r, t) cos 9, sin y(r, t) sin 9, cos y(r, t)). (3-4) 

Assume dm) is true. Then a straightforward calculation as in Section 2 implies that ip solves 
(j2.12p . To see (13.4|) . first note that we can always write 

d(x, y,t) = (^ sin ip{r, 9, t) cos x(r, 9, t), sin y(r, 9, t) sin x(r, 9, t), cos ip{r, 9, t)) (3-5) 

for a pair of functions ip and x- For a G [0, 27r], let Ra G SO(3) be the rotation map in of angle 
a with respect to the 2 -axis, and Sa G SO(2) be the rotation map in of angle a with respect to 
0. Since the initial and boundary values do of d are axisymmetric, i.e., for any a G [0, 27r], 

do{Sa{x, y)) = Rao do(x, y), (x, y) G Bf, 

we can check that R~^ o d[Sa{x,y),t), Va G [0,27r], also solves the initial and boundary value 
problem of the drifted harmonic map equation (|2.12l) . By the uniqueness of solutions to (|2.12l) . 
this implies that 

d{x, y, t) = R~^ o d(5«(x, y),t) , Va G [0, 27r], (x, y) G Bf, 

or equivalently, 

d(x, y,t) = (^ sin ip{r, t) cos{r]{r, t) + 9), siny(r, t) sin(y(r, t) + 9), cos ip{r, t)) (3-6) 
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for some functions (p{r,t) and rj(r,t). Direct calculations, similar to Section 2, imply that ip and rj 
satisfy 


(ft + ripr = {<Prr + —) “ + \) sin <P cos if, 

sin ip{rit + rr]r) = sinip{^r]rr H-) + 2(prr]r cos ip. 


(3.7) 


Now we need to show 


Claim 1. r]{r,t) = 0 for {r,t) G [0,1] x [0,to)- 


Note that when sin(^(r, t) = 0, we can simply define r]{r,t) = 0 because any value of r]{r,t) gives 
the same value of d by (I3.6p . Set 


1/ = I (r, t) G [0,1] X [0, to) I sin ip{r, t) / o|, 

and dehne 

h{r,t) = ri{r,t) sin ip{r,t), {r,t)&U. 

Then by (I3.7I) . h{r,t) satisfies 

ht + rhr = hrr + ^+ [ipl - {r]l + cos^ ip\ h, 


under the initial and boundary conditions; 


(3.8) 


h{r,t) = 0, on dU n {0 < t < to}- (3.9) 

Since ipr is bounded in U, replacing h by he~^ for a sufficiently large c > 0 we may assume the 
coefficient of h in (j3.8p is negative. Therefore by the maximum principle, we conclude that h = 0 
and consequently t] = 0 on U (see [2] Lemma 2.2, or [12] and I37|). Hence Claim 1 holds and the 
proof of Theorem 3.1 is complete. □ 


4 Existence of global smooth solutions and proof of Theorem 12.5 

This section is devoted to the proof of Theorem 12.51 on the existence of global smooth solutions 
when the initial data cfo satisfies (/>o(0) = 0 and |<))o(^)l < for all r G [0,1]. The proof is motivated 

by m- 

We start with 

Lemma 4.1 (maximum principle) For any smooth solution ip{r,t) : [0,1] x [0,T) —?• R to the 
drifted harmonic map equation (I2.12p . if 

(/?o(0) =0, -TT < ipo{r) < vr, V r G [0,1], 


then 

— vr < ip{r, t) < TT, V r G (0,1) and t G (0, T). (4.1) 

Proof. We only present the proof for ip < tt, since the other half of the inequality (4.1) can be 
proved similarly. Set ip = tt — ip. Then ip satishes 

ipt + ripr= Prr-\ - \-pi[r,t)(p, (4.2) 

r 
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where 


and 


Pi{r,t) 


sin(2(^) 
2(pr‘^ ’ 


(f{r, 0) > 0, ip(l, t) > 0 and (^(0, t) = tt. 


Thus for any t G (0, T), there exists a small ri G (0,1) such that (p{r,T) > 0 for any (r, r) G 
(0,ri] X (0,t]. On (ri,!) x (0, t], by the fact that pi(n,r) is bounded, we can prove (^ > 0 by the 
standard maximum principle (see [l2] or [37]). Therefore, we conclude that (^ > 0 or equivalently 
(/? < vr on (0,1) X (0, T). □ 


Lemma 4.2 (comparison principle) Suppose the functions p, f and g are smooth solution, subso¬ 
lution and supersolution to (|2.14l) on [0,1] x [0, T) respectively, and 

f{r, t) < ip{r, t) < g{r, t) on ([0,1] x {0}) U ({0,1} x (0, T)). 


Then we have 

f{r,t) < p{r,t) <g{r,t), V (r,t) G [0,1] x [0,r). (4.3) 


Proof. Set f = f — T- Then / satisfies 


ft + rfr > 


frr + — +P2{r, t)f, 


(4.4) 


where 

sin(2/(r,t)) - sin(2(y9(r,t)) 
2r2(/(r,t)-^(r,t)) ’ 

and 

/(r, t) < 0 on ([0,1] X {0}) U ({0,1} x (0, T)). 


For any t G (0,T), there exists a small r 2 G (0,1) such that P 2 ir,T) < 0 on (r, r) G (0, r 2 ) x (0,t). 
Combining with the fact P 2 ir, r) is bounded on (r 2 ,1) x (0, t), we conclude that P 2 ir, r) is bounded 
from above on (r, r) G (0,1) x (0,t). By the standard maximum principle (see [l2] or [37]), we 
conclude that f <0 or p > f on [0,1] x (0, T). Similarly, one can prove 4> < g. □ 


Proof of Theorem 12.51 To prove the existence of global smooth solutions of (I2.12h . we need to 
construct suitable supersolutions and subsolutions to (I2.14p . Denote 


ip{r, c) = 2 arctan 



and 


p{r, c) = 2 arctan 



(4.5) 


for some positive constant c. It is easy to see that p{r,c) and p{r,c) are smooth functions in [0, Ij. 
Direct calculations give 


_ _ _ Pr sin(2(/?) 

ipt T Trr rt 2 

r Ir 


2rc 

g2 _j_ J.2 


>0, rG [0,1]. 


(4.6) 


Thus is a supersolution of (I2.14p . Denote 


rjo = 7r- rnax |(^o(r)l. 
0<r<l 
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By the assumption on (fQ, we have % > 0 and |(/?o| < — 7?o- Since (^^(0, c) = | and <^(0, c) = 

(^o(O) = 0, we can find a sufficiently small c > 0 such that 

Tp{r,c) >^pQ{r) (4.7) 

for any r € [0,1], with equality iff r = 0. Similarly, we can prove that ip{r,c) is a subsolution to 
(j2.14p and 

ipo{r) > ^{r,c) m[0,l] (4.8) 

for a sufficiently small c > 0. By lemma [121 we can conclude that Tp[r,c) > (p{r,t) > ip{r,c) for 
r E [0,1] and t > 0. 

Suppose T > 0 is the maximum time interval for cp. For any ro E (0,1) and ti E (0,T), by the 
standard regularity theory of parabolic equations, we can prove 

ll‘'^llc''([ro,i]x(ti,T)) - ^^0 > VA: > 1. (4.9) 

This implies that the possible singularity of the solution cj) can only happen at r = 0. Suppose 
T < +00. Then 

d{x, y, t) = (sin (p cos 6, sin (p sin 9, cos ip) 

blows up at ((0,0),T). By the standard blowing up argument (cf. [H] Theorem 1.3), there 
exist {xm,ym) (0,0), tm \ T, r-m i 0 as m ^ Too and a nonconstant smooth harmonic map 
w : ^ such that 


dm{x,y,t) := d{rmX,rmy,tm + rl,t) ^ uj in x (-oo, 0]). (4.10) 

Since, by lemma 021 

-TT < ip{r,t) < TT, r E [0,1], tm<t<T, 
we have that the third component, of w satisfies 

— 1 < oj^{x,y) < 1 for any {x,y) E 

This implies that u E C'°°(M^,§^) is a nontrivial harmonic map, with finite energy and degree zero, 
which is impossible. Therefore, we conclude that T = +oo. This proves Theorem 12.51 □ 


5 Finite time singularity and proof of Theorem 12.6 

In this section, we will modify the construction by [3] to show the existence of solutions with finite 
time singularity of (|2.12l) for some suitably chosen initial data ipo with (/?o(0) = 0 and |¥5o(l)l > 


Proof of Theorem 12.61 We adopt some ideas of [5] to construct suitable barrier functions from 
below. Without loss of generality, we assume that </?o(l) > tt- Denote the scalar-valued tension 
field by 


r((^) := 


^rr T 

r 



There are a family of smooth solutions of t((^) = 0, with 0(0) = 0 and 0(r) > 0, given by 


0(r, /3) = arccos 


//32-r2\ 


or 


0(r, /?) = 2 arctan 



(5.1) 


for /3 > 0. 
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In order to handle the drift term r(pr in the equation (12.1411 . we need to reparameterize (3 by 


e'/3 in (15.1|) : 


4>{r, j3,t) = arccos 


g2t/j2 _ ^2 
^2tp2 J.2 


or 4){r, I3,t) = 2 arctan 




tR ' 


(5.2) 


which is also a solution of T{(j)) = 0 for any t G [0, +oo). For any e G (0,1) and /r > 0, let a = 1 + e 
and 

/ J.0 

6{r, n,t) = 2 arctan 




Then satisfies 


I 0 


1 „ oP sin 9 cos 6 


= 0 . 


r 

Choosing /i large enough so that 9{r,fi,t) is small enough and 

cos 9{r,fj,,t) > 

for any r G [0,1] and t G [0, +oo). We will look for a subsolution / of (|2.14p in the form: 

/(r, t) = (j){r, + 9{r, /i, t), 

where /3(t) solves 

1 f = -fc-”;?'. 

\/9(0) = A, 

where S and Po are positive constants to be determined later. 

Claim 2. If 0 < 2/3q ^ < (5(1 — e), then there exists 0 < Tq < +oo such that 

lim /r(0, t) = +00. 
ttTo 

To see (|5.8I) . we solve the ordinary differential equation (|5.7p and obtain 

ll-e _ ~ /„-2t I ol-e 




-{e-^^-l) + P^ 


Set 

Then 

By (15.1211 below, we have 

so that 


Tn := ^ In 


<5(1 - e) - 2/3, 


l — £ 


> 0 . 


P{t) —0 as t ^ Tq 


„ , . 2/3e* 2a/rr^e“* 

Jr{r,t) = 0.^0 . —^ + 


friO,t) = 


2 

+ 00 , as t —)• Tq . 


e^P{t) 


Claim 3. There exists 5 > 0 such that f is a subsolution of 112. 12\) . i.e., f satisfies 

ft + rfr < r(/). 


(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


(5.8) 


(5.9) 
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In fact, by (15.11) and (15.41) . we have 


r(/) =— [ sin (p cos p — sin((/) + 9) cos{(f) + 9) + sin 9 cos 0] 
[a^ sin 9 cos 9 — cos( 2 (/i + 9) sin 0] 

1 


> — [(1 + e) sin 6 ( — cos(2(/) + 9) sin 9] 
esm9 

— 


e 2|Ue 


—at^a 




> 


2fie 


—at^ 


1 


- ;U2 + 1 

where we have used (j5.5j) and (j5.3p . From (15.Ij) and (15.3j) . we have 

2 re*(/? + j3t) 2 a^r“e“* 


and 


ft = - 


fr = 


g2t/32 ^2 

2/3e* 


_j_ 


2a 


+ 


2 a/rr' 


a—1 


^2tp2 _|_ ^2 ij2^2at _j_ ^2a 

Combining (15.111) with (15.121) and using (15.7p . we obtain 

2re^j3t 25re~*/3^ 


ft + rfr = -- 


g2t^2 _j_ ^2 e 2*/32 + r2 ' 

To prove Claim 3, it suffices to verify the following inequality: 


25re ^ 27 /re “*e 


< 


—1 


Let s = 


g 2 t /52 _|_ J.2 — /i2 

Then ()5.14l) is equivalent to 

^2-e 


< 


7 /re 


1 + ^2 (5(/r2 + 1) 

^2-e 


, Vs > 0. 


(5.10) 


(5.11) 

(5.12) 

(5.13) 

(5.14) 


(5.15) 


It is easy to check that the function 
if we choose 


1 + s2 


5 < 


has a maximum M (e) depending only on e. Therefore, 
7 /re 


M(e)(/r2 + l)’ 
then ()5.14p holds and hence the Claim 3 follows. 

Claim 4. For sufficiently large /r > 0, there exists (po G C'°°([0,1]), with VJo(O) = 0 (/5o(l) > tt, 

such that 

f{r, t) < (po{r) on ([0,1] x {0}) U ({0,1} x (0, Tq)). (5.16) 

Since v^o(l) > tt, we can choose a sufficiently large /r such that 0(1,/r,t) < </^o(l) — tt for any t > 0. 
This, combined with 0 < (j){l,l3,t) < vr, implies 


/(I, t) = 4>{l, (3, t) + 0(1, /r, t) < (^o(l), VO < t < To. 
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It is clear that we can find some initial data </?o G C'°°([0,1]), with <^o(0) = 0 and |<?io(l)| > tT) such 
that ipo{r) > /(r,0) for any r G [0,1]. Hence (I5.16p holds. 

It follows from Claim 3, Claim 4, and lemma 021 that 

/<(/?, in [0, 1] X [0,To). 

This, combined with /(0,t) = = 0, implies that (pr{0,t) > fr{0,t) for 0 < t < Tq. Hence, by 

Claim 2, we have 

t) —>■+ 00 , as t—)-T q". 

This completes the proof of Theorem 12.61 □ 


Remark 5.1 It is clear that the Claim 3 doesn’t hold if the time relexation constant 7 = 0. In 
fact, it is a very challenging question how to construct solutions with finite time singularity to 
(j2.6j) when 7 = 0 . Such an example would shed lights to the question whether the equation of 
viscoelastic flows admits finite time singularity in dimension three. 


We finish this section with two remarks on Theorem 11.1 


Remark 5.2 Under the same assumptions as Theorem ll.il does the stress tensor S := Vd© Vd — 
-jVdpIa blow up as t approaches Tq? 

By calculating the stress tensor for the corresponding subsolution, we conjecture that S doesn’t 
blow up. Here we sketch the calculation. Assume {u,d,P) is given by Theorem ll.il and / is given 
by Section 5 above. Set 

d = (sin / cos 6, sin / sin 9, cos /). 

Then direct calculations imply 

^ cos ^ sin^ 9 (/^- ^ ) cos 9 sin 9 ^ 


Vd 0 Vd = 


( f 2 f \ a ■ a f‘2 ■ 0 “^ I f 2 a 

y [Jr -2—j cosdSind smd H-^—cos d 


and 


At r = 0: since 


|Vd|2 = /2 + 


sin^ / 


sin/ 


= fr = 


e , it follows that 


m 

Vd 0 Vd 


1 

f /301) 

0 

\ 

r=0 


4 



WW) 

/ 


.-21 


and 


|Vdp 


.-21 


r=0 


Therefore the stress tensor for d at r = 0 is 


(5.17) 


(5.18) 


(5.19) 


(5.20) 


5 := (Vd© Vd - -iVdpIa) 


r=0 


0 0 
0 0 


(5.21) 


where I 2 is the identity matrix of order 2 . 
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Remark 5.3 Under the same assumptions as Theorem 11.11 does the solution {u,d,P) satisfy the 
energy dissipation inequality 

iiil f ^ ^ ^ |Ad+ iVdpdp) < 0, (5.22) 

for 0 < t < To? 

Since u{x,y,z,t) = uo{x,y,z) = {x,y,—2z) in Q, = Bf x [0,1], direct calculations imply 


-/ 

dt Jn 


up = 0, / Au • u = 0, 


and 


f (n • Vn, n) = [ u- V(^) = f (x^ + - 82:^) = - 

J Q, J Q. " J Q. 

/ jVup = Ovr. 


13 


-vr, 


O) 1 can be rewritten as 


v(y + ^ + 222+p) =-(Ad,Vd). 


Multiplying (j5.23l) by u and integrating over U yields 

—TT — f Pu ■ u = I {Ad,u -Vd). 

6 Jan Jn 

Since the boundary condition of d on is given by 

dd 

d = do on dB\ x [0,1], ^ = 0 on B^ x {0,1}, 
multiplying (ll.ip -i by Ad and integrating over fl we would have 


t/ t/ 

Adding (5.24) and (5.25) together, we have 


2j|2 


in 


|Ad+ |Vdpd| 


Jtili + L + |Ad+ jVdpdp 


13 f 

= Ott H-TT — / Pu ■ V. 


Therefore (I5.22p holds iff 


I an 


13 


6 


/ 

Jan 


Gvr H-TT — / Pu ■ u < 0. 


However, it is unclear whether (|5.27p holds for the solution by Theorem 11.1 


(5.23) 


(5.24) 


(5.25) 


(5.26) 


(5.27) 
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6 Finite time singularity for generic initial data and proof of The¬ 
orem 11.31 

This section is devoted to another construction of finite time singularity of (II.ip for more generic 
initial data, in which the solution satisfies the energy dissipation inequality (|5.22|] . 

First we recall the following result on the existence of local smooth solution to (jl.ip and (jl.Sp . 
whose proof can be found in |26j . 

Lemma 6.1 For {uo,do) E M^) x there exist Tq = ro(rio,fio) > 0 and a 

unique smooth solution {u, d) E x [0,To),M^ x S^) to the system m along with the initial¬ 

houndary condition (Egi. Moreover, the energy dissipation inequality 115.22\) holds for 0 < t < Tq. 

Now we would like to proceed with the proof of Theorem 11.31 as follows. 

Proof. Assume Tq > 0 is the maximal time interval for the short time smooth solution (u, d) by 
lemma [6Tl We want to show 

Claim 5. If eo > 0 is sufficiently small, then Tq < 1. 

We argue by contradiction. Suppose that Claim 5 were false. Then for any e > 0 we can find 
{uo,do) E X Cf^{Bf,§‘^) such that 

(a) do is not homotopic to e relative to dBf, 

(h) E{uo,do) < _ 

and a smooth solution {u,d) E C°°[Bf x [0,1],M^ x S^) to (11.11) and (II.8p . Integrating (15.221) over 
t yields that {u, d) satisfies the energy inequality: 

E{u{t), d{t)) + f [ (|Vu|2 + I Ad + I Vdpdp) < E{uo,do) < (6.1) 

Jo JbI 

for all 0 < t < 1. Applying Fubini’s theorem to (|6.1I) . we find that there exists ti E (^, 1) such that 

E{u{h), d{ti)) + [ (|Vu(ti)|2 + |Ad(ti) + |Vd(ti)|2d(ti)|2) < 8e2. (6.2) 

JbI 

From (16.21) and e-apriori estimate ()6.6I) of Theorem 16.21 below, we conclude that there exists a 
universal C > 0 such that 

[d(ti)]^i^^3^ < (6.3) 

Thus d{ti){B\) C B^^{e) n and hence d{ti) is homotopic to e relative to dBf, provided e > 0 

is chosen to be sufficiently small. Since d E C°°[Bf x [0,ti],S^) and d = e on dBf X [0,ti], we see 
that d{ti) is homotopic to do relative to dBf and hence do is homotopic to e relative to dBf. This 
contradicts the assumption (a). Hence Claim 5 is true. This completes the proof of Theorem 1 1.31 □ 

Now we need an e-apriori estimate on approximate harmonic maps from Bf to which can 
be proved by suitable modifications of the arguments by Ding-Wang [9] and Lin-Wang [29] . 

Theorem 6.2 There exist e > 0 and C > 0 such that if d G Cf^{Bf,Ef) satisfies 

E{d) :=l [ |Vd|2 < e2, (6.4) 

^ JBf 
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and 


Then 



Ad+ iVdpdp 


< 


(6.5) 




< CVi. 


( 6 . 6 ) 


Proof. The proof is based on suitable modifications of that by Ding-Wang [9] and Lin-Wang [29]. 
For the completeness, we sketch it here. We divide the proof of estimate ()6.6p into two lemmas. 


Lemma 6.3 Under the same assumptions as Theorem, \6.i^ for any fixed Sq G (e, 1) we have that 

[dl 1 , .<CV^, (6.7) 

holds for any xq G and 0 < ro < do- 

Proof of Lemma 16.31 In order to show (j6.7p . set the tension field of d by 

T(d) = Ad + iVdl^d ; Bf E^. 

We will first establish a modified energy monotonicity inequality for approximate harmonic maps 
in dimension three. From (16.5p . we see that tension field r(d) G L‘^{Bf), and 




L2(S3) ^ «• 


( 6 . 8 ) 


Now we have 

Claim 6. For any xq G Bf and 0<r<i?<l — |xo|, it holds 


-7. 


Br{xo) 


iVdp < 8R 


-1 


IBii{xo) 


\Vd\‘^ + 8R 


L 


Br{xo) 


|r(d)p 


(6.9) 


To see (j6.9h . we assume for simplicity that xq = 0. Multiplying the approximate harmonic map 
equation 

Ad+|Vd|2d = r(d), in Bf, (6.10) 

by X-Vd, integrating the resulting equation over for 0 < r < 1, and applying the same argument 
as |29| lemma 5.3, we obtain 



(T(d),x-Vd))) = r ^ f 1^1^ 
^ JdBrio) dr 



( 6 . 11 ) 


Thus (|6.9p follows by integrating (|6.1ip over r and Holder’s inequality. 
It follows from (16.4p . (16.5p . and (16.9h that 


sup 


{r ^ f |Vd|^ : X G Bf_s^, 0 < r < do} 


' Br(x) 


< 8d(7^ / |Vdp + 8do [ |T(d)|2 < Ce^{l + 5n^) < Ce. 
IBf JBl 


( 6 . 12 ) 


Next we have 
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Claim 7. For any a E (0, ^], there exists 6 q E (0, such that for any xq E Bfso 0 < r < (5o 
it holds 


{dor 




fBr{xo) 


\r{d)\\ 


|vd|2 <f |vd|2 + 

J Br(xo) 

To show (|6.13l) . first recall the Morrey space for U C 

m2-2([/) := // e lL(C/) I \\f\\U2(^u) = sup r-i [ \f\^ < +oo|. 

Br(x)CU JBrix) ’ 


(6.13) 


From ()6.12p . we have that Vd E M^’‘^{Bsq{x)) for any x E and 


llv^ll 


M^,HBs^ix)) 


< Ce, Vx E Bf_ 


(^0* 


(6.14) 


For Xq E Bf_g^ and 0 < r < 1 — (5o, let d E Fr^(M^,M^) be an extension of d such that d = d in 
Br{xQ), |d(x)| < 2 for X E and 


f II^'^IL 2 (K 3 ) < ^\\^A\L‘i{Br{xo)Y 
III^'^IIm 2.2(IR3) < ^||^^|lM2,2(B^(a;o))- 

By the Helmholtz decomposition, there exist G E dd^(M^) and H E L^(M^,]R^) such that 
j Vd X d = VG + dd, V • dd = 0 in R^, 

ll|VG|| L2(]R3) + II'^IL2 (r3) < g||v<!||„,.., < c||vd|| L'2{Br{xo))- 


(6.15) 


(6.16) 


Applying the Poincare inequality, the duality between the Hardy space dd^(R^) and BMO(R^) (see 
m, [15], or [29]), we can estimate 


IddP = 


H-iVdxd- VG) = / H-Vdxd 


- ^11-^ ■ ^'^II«1(K3)II'^IIbMO(R3) 

- ^l|■^llL2(R3)||^^IL2(R3)||Vc^||M2,2( 

- ^\\"^A\L‘2{Brixo))\\'^A\M‘2,2(Br{xo)) “ ^ A\L'^ (Br (xo))' 

From ()6.10p . G solves 

AG = T(d) X d in ddr(xo). 

By the standard L^-estimate, we have that for any 0 < 0 < 1, 


f \VG\'^<Ce^f iVdl^ + r^f |T(d)|2. 

J Bgr{xo) J Br(XQ) J Br{xo) 


I Bgr{xo) 

Putting p6.17p and (I6.19P together, we obtain 


|Vdp < G(02 + 0-ie)i /" |Vd|2 + G0-V 

dT' JBgr{xo) JBr(xf)) 


JBr(xo) 


' Bq- (^XQ ) 


|T(d)P 


(6.17) 

(6.18) 

(6.19) 

( 6 . 20 ) 


for any xq E Bf_g , 0 < r < 5o, and 0 < 0 < 1. 
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For any 0 < a < ^, first choosing 6 = 9o ^ (Oj 1) such that C9q < ^9q^ and then choosing 


e = Co > 0 such that C9q < ^9^°^, we obtain (16.131) . It is standard that iterations of (I6.13|) imply 
that 

1 /■ „ r ^ i2 


1 


® JBs{xo) 


|Vd|2 < (-)2« /" iVdp + Cs/" 
^ JBrixo) Jb 


' Br (^0 ) 


\r{d)\^ 


( 6 . 21 ) 


holds for any xq G and 0 < s < r < dg. It is clear that (|6.7I) follows from (|6.21l) and Morrey’s 

decay lemma (cf. [36]). □ 

Lemma 6.4 Under the same assumptions as in Theorem Ih.M there exists 0 < dg < ^ depending 
only on dBf such that for any xq G dBf and 0 < rg < 2(5g, it holds 




( 6 . 22 ) 


where Bf~{xo) := Br{xo) n Bf for r > 0. 


Proof of Lemma 16.41 The strategy to show (I6.22p is similar to lemma 16.31 that is to establish a 
modified boundary monotonicity inequality and a boundary energy decay property. More precisely, 
we first need 


Claim 8. There exists 0 < dg < depending only on dBf, such that for any xq G dBf and 
0 < r < R < 5o, it holds 


\Vd\‘^<8R-^[ lVd|2 + 16i? 
JBf{xo) J Bf(xo) 


Ib+{xo) 


(6.23) 


For r > 0, set B+ = Br{0) n = {xi,X 2 ,xs) G : xg > O}, T^ = B+ n {x = (xi,X 2 ,X 3 ) G 
: xg = O}, and Sf' = 9B+ \ T^. By the standard boundary flatten argument, there exists 
(5g > 0 depending only on dB\ such that for any xg G dB^ there exists a smooth diffeomorphism 
4>g : B|^ -S^(3;o) such that 4'g(0) = xg, ^oiTso) = n dBf, and 

|^o(a;) - xg| < C\xf, |V4'g(x) - Id| < C|x|, |V2^>g(x)| < C, 'dx £ B^^, (6.24) 

where Id denotes the identity map in 

Set g = 4’g((7g) in B^, pull-back of the euclidean metric go = dx^ in i?^(xg). Consider d = 
do <I>g : B^ —>■ Then it is not hard to check that d : (B^,gf) —>■ is an approximated harmonic 
map, i.e., 

Agd-b IVd|^d = T(d ) in B((; d = e on T^^, (6.25) 

and 

|r(d )|(x) < C'|r(d)|(4'g(x)), Vx G B^^, (6.26) 

where is the Laplace operator with respect to g, and |Vd|g is the Dirichlet energy density of d 
with respect to g. Note that by (|6.24l) . g satisfies 

\g{x) - go\ < C\x\^ |V 5 (x)| < C|x|, \V^g{x)\ < C, Vx G B+. (6.27) 

Based on (I6.27p. (|6.25p. and (|6.26l). we mav for simolicitv assume that xn = 0 and d : (“1 , do) ^ 
is an approximated harmonic map with tension field r(d). Since d = e on it is easy to see 
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X • Vd = 0 on Tsq. For 0 < r < 5o) multiplying (|6.10p by x -Vd and integrating the resulting 
equation over B+, we get 


1 

2 



r 

2 




This implies 



Integrating (16.28p over 0 < r < i? < do and applying Holder’s inequality, we obtain 


(6.28) 


r-^ [ |Vd|2 < 8i?"^ / |Vd|2 + 16i? / |T(d)p. 

JB+ ^8+ Jb+ 

This gives p6.2,8p . 

Next we need 

Claim 9. For any a G (0, ^), there exists 6 q G (0, 5 ) sueh that for any xq G dB^ and 0 < r < 60 , 
it holds 

{0or)~^ [ |Vd|2 < / |Vdp + Cv/ |T(d)p. (6.29) 

JBf(xo) JEfixo) 

For simplicity, we again assume xq = 0 and H^(xo) = B^. The proof is similar to [9] lemma 3.3, 
and we only sketch it. To obtain (I6.29p . we perform suitable extensions of d to B^^ as follows. Let 
d : B^p be the extension of d that is even with respect to X 3 , and let r : B^p —>■ be the extension 
of r(d) that is odd with respect to X 3 . Define w : B^^ —by 


w{x) = 


{d — e)(x) if X = (xi, X2,X3) G B^p and X3 > 0, 

-(d-e)(xi,X2,—X3) if X = (xi,X2,X3) G B^^ and X3 < 0. 


(6.30) 


For 0 < r < do, let r/ G Cj 


be even with respect to X 3 , 0<t/< 1, r/ = lin Bn, r] = 0 


outside Br, and \Vr]\ < Cr To proceed with the proof, we need 

V • (Vw X d ) = T X d in D'(B^(,). 


(6.31) 


To see (j6.3ip . let (p G C^(B5p) and write 4> = 4>e + 4>o, here 4>e i4>o) is even (odd) with respect to 
X 3 respectively. Then (|6.3ip follows from 

/ Vwxd-Vf = 2 Vwxd-V4>o = —2 / V ■ {Vd x d) ■ fo 

Jbs„ JBt JBf 


= -2 


j T{d) 

JBt 


I 


X d ■ 4>o = — / T X d ■ 4>. 


It follows from (j6.9ll and (j6.23l) that Vrc G M‘^’‘^(Msg), and 


(6.32) 


From p6.3ip . it is easy to check that V • {rfVw x d ) G L^(M^), and 


(6.33) 
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Observe that 


\Vd\‘^<2[ X dl"^ = [ r]‘^\Vwxdf 

JkI Jrs 


= — / V{rj‘^Vw xd)-{wxd) + j 7f[{yw x d) x Vd — A] • w + A J jfw 


= / + // + III, 


where 


A = 


rfiVw X d) X Vd 
/k3 ri^ 


By the Poincare inequality, we have 

|m| < |A| / \w\<C[r-^ j \Vd\^Y j |Vdp 
^ 8 + JB+ Jb+ 

Applying (j6.3,'ljl . I can be estimated by 

III < V • irt^Vw X d) w 

^ ' L2(R3) L2{B+) 

d — e 


— II^'^IIi, 2 (b+) + ||T((i)||^ 2 (B+) 


< 7 

for any 7 G (0,1). Let w : 


I 2 || / tmi 2 

_j_ ^ 


L2(B+) 

+.)+C' 7 -V- 2 / |d-e| 
’■ '2 7™+ 


—>■ be an extension of w such that 

||Vui||^2,2(R3) < ^||^^|Im 2 . 2 (B 5 ^)- 

Then, from (j6.32|] and the Poincare inequality, we have 

MbMO(R3) — ^II^^IIm2,2(K3) < ^II^^IIm2.2(B^^) — C\/e. 

We use the duality between T-L^{M.^) and BMO(M^), similar to [S] lemma 3.3 and lemma 
estimate II hy 


\II\ < C 
< C 


rj [(Vw X d) X Vd A] [H bmo(R3) 


II^^IL2(b+) + ^ 11^ ■ ^ '^)|Il2(b+) II^^IIm2.2(bj^) 

< C^/e ||Vd||^2(B+) + ^ II'''(^)|Il2(b+) • 

Putting these estimates together, we obtain 

1 


^ / |Vd|^ < 

^/2 7b+, 

r/2 


C{V~e + l)\- [ \Vdf + r[ \r{d)f]+Cy-^\ f \d 
L?' 7b+ 7b+ -I 7b+ 


— e 


To deduce (16.290 from (j6.35p . we need 

Claim 10. There exists Oq = 9o{e) such that for 0 < r < 60 , 


{9or)^ Ie+ 

Oor 


-e |2 < CB'^- 


7b+ 


|Vdp 


(6.34) 
2 . 6 , to 


(6.35) 

(6.36) 
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We prove (I6.36P by contradiction. Suppose that it were false. Then for any 6 E (0,1), there exist 
Cfc —)• 0, 0 < Tfc < (5o, efc E and dk E such that dk = on 


„-i 


but 


|V4|^ = el and f |r(4)|^ < Cfc, 

3 ' K 

[ \dk-ek\^>ke^- [ |V4I". 

Jb+ Jb+ 


(6.37) 


(6.38) 


Define 4(3;) = ^ —>■ S^. Then we have 




and 

but 


4|-r=0, / |V4p = i,, 

A4 = fk ■= -efc|V4p4 + e^V(4) 0 in T^(B]^), 

i [ \dkf > kB'^. 

JMt 


(6.39) 

(6.40) 

(6.41) 


From (j6.39p . {4} C 77^(B+,M^) is bounded. We may assume that dk ^ d weakly in i7^(B^ 
strongly in (B^). It follows from (|6.39p and (16.4011 that 


d\T^=0, /jVdp 

4b+ 


< 1, 


and 


Ad = 0 in 


"1 ■ 


By the standard theory of harmonic functions, we have that for any 6 E (0,1), 

1 [ \df < CB\ 

JMj 

Since dk ^ d in L^(B)^), (I6.44p contradicts to (|6.41l) . Hence Claim 10 is proven. 
Putting (I6.36P into ()6.35p . we can obtain that 


(6.42) 

(6.43) 

(6.44) 


1 

Bor 


[ |Vd|2 < C{Ve + j)\^ [ |Vd|2 + 4r 


k(d)p 


200 r 
2 


"’200’- 


JB+ 

< C[(v^ + 7)C+7"'0o] 4 /jVd|2 + C4r / jT(d)p 

r Jb+ Jmj 


< CB, 


2a 


4/ |Vd|2 + C4r/ \T{d)\\ 
jvt 4B+ 


(6.45) 


for any 0 < a < i, provided we choose 7 = dg ^ — G‘o^^°‘■ Hence (I6.29P is proven. 
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(6.46) 


It is standard that iterations of (j6.29l) imply that for any a E (0, i 


2n 


-[ 

S JB+ixo) 


|Vd|2 < / |Vdp + Cs/ |r(d)| 

r Jb+{xo) Jb+{xo) 


holds for any xq E dBf and 0 < s < r < (5o. Combining (16.461) with (16.211) . we can obtain that for 
any 0 < a < ^, 


1 

s 


I Ba{xo)nBf 


|Vd|2 < / |Vdp + Cs / 

^ JBr{xo)nBf Je 


I Br(xo)nBf 




(6.47) 


holds for any xq E Bf and 0 < s < r < (5o- 
implies that for any a E (0, ^), 


This, combined with Morrey’s decay lemma (see H) 




C' 


(i?3) 


< C 


l|Vd|| 


L\Bl) 


+ |k(c?)|| 


L2(Bf) 


(6.48) 


It is well known that we can use the equation (|6.25p to improve the estimate (16.481) to the case that 
a = ^. Hence lemma [63] is proven. This completes the proof of Theorem 6.2. Thus Theorem 1.3 
is proven. □ 


Acknowledgements. The first author is partially supported by NSF grants DMS 1412005 and 
DMS 1159937. The second author is partially supported by NSF grants DMS 1065964 and DMS 
1159313. The third author is partially supported by NSF grants DMS 1412005, DMS 1216938 and 
DMS 1159937. The fourth author is partially supported by NSF grant DMS 1522869 and NSFC 
grant 11128102. 


References 

[1] K. C. Chang, Heat flow and boundary value problem for harmonic maps. Ann. Inst. H. Poincare 
Anal. Non Lineaire, 6 (5) (1989), 363-395. 

[2] K. C. Chang, W. Y. Ding, A result on the global existence for heat flows of harmonic mpas 
from into S'^. Nemantics, J.-M. Coron et al. ed., Kluwer Academic Publishers, 1990, 37-48. 

[3] K. C. Chang, W. Y. Ding, R. Ye, Finite-time blow-up of the heat flow of harmonic maps from 
surfaces. J. Diff. Geom., 36 (1992), 507-515. 

[4] Y. M. Chen, M. Struwe, Existence and partial regularity results for the heat flow for harmonic 
maps. Math. Z., 201 (1) (1989), 83-103. 

[5] Y. M. Chen, F. H. Lin, Evolution of harmonic maps with Dirichlet boundary conditions. Comm. 
Anal. Geom. 1 (1993), no. 3-4, 327-346. 

[6] Y. M. Chen, W. Y. Ding, Blow-up and global existence for heat flows of harmonic maps. Invent. 
Math. 99 (1990), no. 3, 567-578. 

[7] J. M. Coron, J. M. Ghidaglia, Explosion en temps fini pour le flot des applications harmoniques. 
G. R. Acad. Sci. Paris, 308 (1989), 339-344. 


25 














[8] H. J. Dong, Z. Lei, On a family of exact solutions to the incompressible liquid crystals in two 
dimensions. arXiv:1205.3697, 2012. 

[9] S. J. Ding, C. Y. Wang, Finite time singularity of the Landau-Lifshitz-Gilbert equation. Int. 
Math. Res. Not. IMRN 2007, no. 4, Art. ID rnm012, 25 pp. 

[10] J. L. Ericksen, Hydrostatic theory of liquid crystal. Arch. Ration. Mech. Anal., 9 (1962), 371- 
378. 

[11] L. C. Evans, Partial regularity for stationary harmonic maps into spheres. Arch. Ration. Mech. 
Anal. 116 (1991), 101-113. 

[12] A. Eriedman, Partial differential equations of parabolic type. Prentice-Hall, Englewood Cliffs, 
New Jersey, 1964. 

[13] J. Grotowaski, Harmonic map heat flow for axially symmetirc data. Mann. Math., 73 (1991), 
207-228. 

[14] J. Grotowaski, Finite time blow-up for the harmonic map heat flow. Calc. Var. Partial Differ¬ 
ential Equations, 1 (1993), 231-236. 

[15] F. Helein, Harmonic maps, conservation laws and moving frames. 2nd ed. Cambridge Tracts 
in Mathematics, 150. Cambridge University Press, Cambridge, 2002. 

[16] J. Hineman, C. Y. Wang, Well-posedness of nematic liguid crystal flow in Arch. 

Ration. Mech. Anal., 210 (2013), 177-218. 

[17] M. C. Hong, Global existence of solutions of the simplified Ericksen-Leslie system in dimension 
two. Calc. Var. Partial Differential Equations, 40 (2011), no. 1-2, 15-36. 

[18] M. C. Hong, Z. P. Xin, Global existence of solutions of the liquid crystal flow for the Oseen- 
Frank model in Adv. Math. 231 (2012), no. 3-4, 1364-1400. 

[19] M. C. Hong, J. K. Li, Z. P. Xin, Blow-up criteria of strong solutions to the Ericksen-Leslie 
system in M^. Comm. Partial Differential Equations 39 (2014), no. 7, 1284-1328. 

[20] E. Hopf, Uber die Anfangwertaufgaben fur die hydromischen Grundgleichungen. Math. Nach., 
4 (1951), 213-321. 

[21] T. Huang, C. Y. Wang, Blow up criterion for nematic liquid crystal flows. Comm. Partial 
Differential Equations, 37 (2012), 875-884. 

[22] S. Leonardi, J. Malek, J. Necas, M. Pokorny, On axially symmetric flows in Zeitschrift fiir 
Analysis und ihre Anwendungen, 18 (1999), no. 3, 639-649. 

[23] J. Leray, Sur le mouvement d’un liquide visqueux emplissant I’espace. (French) Acta Math. 63 
(1934), no. 1, 193-248. 

[24] F. M. Leslie, Some constitutive equations for liquid crystals. Arch. Ration. Mech. Anal., 28 
(1968), 265-283. 

[25] F. H. Lin, Nonlinear theory of defects in nematic liquid crystal: phase transition and flow 
phenomena. Comm. Pure Appl. Math., 42 (1989), 789-814. 


26 


[26] F. H. Lin, J. Y. Lin, C. Y. Wang, Liquid crystal flows in two dimensions. Arch. Ration. Mech. 
Anal., 197 (2010) 297-336. 

[27] F. H. Lin, C. Liu, Nonparabolic Dissipative Systems Modeling the Flow of Liquid Crystals. 
Comm. Pure. Appl. Math., Vol. XLVIII, (1995), 501-537. 

[28] F. H. Lin, C. Y. Wang, Global existence of weak solutions of the nematic liquid crystal flow in 
dimensions three. Comm. Pure Appl. Math, to appear. 

[29] F. H. Lin, C. Y. Wang, On the uniqueness of heat flow of harmonic maps and hydrodynamic 
flow of nematic liquid crystals. Chinese Ann. Math., 31B (6) (2010), 921-938. 

[30] F. H. Lin, C. Y. Wang, Recent developments of analysis for hydrodynamic flow of nematic 
liquid crystals. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 372 (2014), no. 
2029, 20130361, 18 pp. 

[31] F. H. Lin, C. Y. Wang, The Analysis of Harmonic Maps and Their Heat Flows. The World 
Scientific, 2008. 

[32] F. H. Lin, P. Zhang, Global small solutions to an MHD-type system: the three-dimensional 
case. Comm. Pure Appl. Math. 67 (2014), no. 4, 531-580. 

[33] F. H. Lin, T. Zhang, Global Small Solutions to a Complex Fluid Model in Three Dimensional. 
Arch. Ration. Mech. Anal., 216 (2015), 905-929. 

[34] F. H. Lin, P. Zhang, L. Xu, Global small solutions to 2-D incompressible MHD system. 
arXiv: 1302.5877, 

[35] P. L. Lions, Mathematical topics in fluid mechanics. Vol 1. Incompressible models. Oxford Lec¬ 
ture Series in Mathematics and its Applications, 3. Oxford Science Publications. The Clarendon 
Press, Oxford University Press, New York, 1996. 

[36] C. B. Morrey, Multiple integrals in the calculus of variations. Die Grundlehren der mathema- 
tischen Wissenschaften, Band 130 Springer-Verlag New York, Inc., New York, 1966. 

[37] M. H. Protter, H. F. Weinberge, Maximum principle in differential equations. Pretice-Hall, 
Englewood Cliffs, New Jersey, 1967. 

[38] M. Struwe, On the evolution of harmonic mappings of Riemannian surfaces. Comm. Math. 
Helv., 60 (1985), 558-581. 

[39] R. Temam, Navier-Stokes equations and nonlinear functional analysis. SIAM, Philadelphia, 
1983. 

[40] C. Y. Wang, Well-posedness for the heat flow of harmonic maps and the liquid crystal flow 
with rough initial data. Arch. Ration. Mech. Anal. 200 (2011), no. 1, 1-19. 

[41] X. Xu, Z. F. Zhang, Global regularity and uniqueness of weak solution for the 2-D liquid crystal 
flows. J. Diff. Eqns., 256 (2012), 1169-1181. 


27 


